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ABSTRACT 

We show that the gas mass fraction of spiral galaxies is strongly correlated with 
luminosity and surface brightness. It is not correlated with linear size. Gas fraction 
varies with luminosity and surface brightness at the same rate, indicating evolution at 
fixed size. 

Dim galaxies are clearly less evolved than bright ones, having consumed only ~ 1/2 
of their gas. This resolves the gas consumption paradox, since there exist many galaxies 
with large gas reservoirs. These gas rich galaxies must have formed the bulk of their 
stellar populations in the last half of a Hubble time. The existence of such immature 
galaxies at 2; = indicates that either galaxy formation is a lengthy or even ongoing 
process, or the onset of significant star formation can be delayed for arbitrary periods 
in tenuous gas disks. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental 
parameters — galaxies: photometry — galaxies: spiral — galaxies: structure 



Accepted for Publication in the Astrophysical Journal 



1 



1. Introduction 



2. The Data 



The evolutionary history of galaxies can be inferred 
in a variety of ways. One approach is to utilize the 
extensive knowledge of the evolution of stars gleaned 
from the HR diagram to construct model galaxies 
(e.g., Tinslcy 1968; Bruzual & Chariot 1993). More 
recently, it has become possible to infer related infor- 
mation in a statistical way from deep redshift surveys 
(e.g., Lilly et al. 1995; Ellis et al. 1996). Piecing the 
various lines of evidence together is not trivial, and 
rather different pictures can emerge from equally valid 
approaches. 

It is essential that the question be well posed. 
What is meant by galaxy evolution? We generally 
imagine that galaxies form out of the gas that emerges 
from the hot big bang. Regardless of the details 
of galaxy formation, the units which we observe as 
galaxies today have converted some fraction of this 
gas into stars. The gas mass fraction fg is thus a nat- 
ural parameter for quantifying the degree of evolution 
of a galaxy. 

This obvious path has not been much pursued for 
the simple reason that most well observed galaxies 
have low gas fractions {fg <^ 0.2). Their evolution 
is well advanced, and the small dynamic range in fg 
provides no useful constraint on models. Indeed, this 
leads to a conundrum in that we appear to live at a 
special time when most disks are near to exhausting 
their gas reservoirs (e.g., Roberts 1963; Kennicutt et 
al. 1994). 

However, those objects which are well studied need 
not be representative. It has recently been demon- 
strated that disk galaxies of low surface brightness are 
quite common (McGaugh et al. 1995a, de Jong 1996, 
McGaugh 1996a; Impey et al. 1996). These galax- 
ies are generally gas rich (Schombert et al. 1992; de 
Blok et al. 1996), and Mui/L increases systematically 
with decreasing surface brightness (Fig. 1). This ra- 
tio depends only on the relative amounts of atomic 
hydrogen and star light, and must be an evolution- 
ary effect. Attempting to understand this evolution 
and its strong dependence on surface brightness is the 
subject of this paper. 

We emphasize that our discussion is restricted to 
disk galaxies, primarily spirals. We do not address 
Elliptical or dwarf Spheroidal galaxies, though we do 
thoroughly examine relevant selection effects. When 
required, distance dependent quantities assume Hq = 
75 kms~^Mpc"^ 



The data consist of optical and 21 cm fluxes. Data 
arc selected to provide as large as possible a database 
for investigating the relation between the gas con- 
tent and optical properties of spiral galaxies, particu- 
larly surface brightness. A prerequisite is that each of 
the quantities of interest are actually measured. Sur- 
face brightness, by which we mean the central surface 
brightness of the disk determined from extrapolation 
to i? = of exponential fits to the observed light 
profile, is not commonly measured and is therefore 
the most limiting factor. Of studies which specifi- 
cally measure surface brightness in the appropriate 
fashion, HI measurements are not always available, 
further limiting the available data. It would be desir- 
able to also have CO measurements as a tracer of the 
molecular gas, but there are essentially no published 
data which meet this additional criterion, and low 
surface brightness galaxies have yet to be detected 
in CO at all (Schombert et al. 1990). Instead, we 
will use the gas phase-morphological type relation of 
Young & Knezek (1989) molecular gas content. A fur- 
ther piece of information required for estimating the 
stellar mass is either an I magnitude or a B — V color 
(see §3.1). In sum, we need to know Mg, fiQ, h, Mhi, 
the morphological type T, and B — V or I. 

The large data sets which meet these requirements 
are those of Romanishin et al. (1982; 1983), McGaugh 
& Bothun (1994), de Jong & van der Kruit (1994), de 
Jong (1995; 1996), and de Blok et al. (1995; 1996). 
Surface photometry has been performed on all galaxy 
images to determine the disk parameters /Uq (central 
surface brightness) and h (scale length), as described 
in the various sources. Most of the data for higher sur- 
face brightnesses (/xo < 23 B mag arcsec"^) are from 
de Jong, whereas most of that for [1q > 23 are from 
our own work. An obvious gap visible in Fig. 1(b) at 
intermediate surface brightnesses between these two 
data sets is filled in Fig. 1 (a) by the data of Roman- 
ishin. Where used, data for the special case of Malin 
1 are taken from Bothun et al. (1987) and Impey & 
Bothun (1989). 

No trends with Mm/L are obvious in any one 
data set alone. Only upon combining them does the 
strong relation in Fig. 1 become apparent. Data over 
this large dynamic range (a factor of 100) in surface 
brightness have not previously been available. 
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2.1. Selection Effects 

The data have been selected to be as inelusive 
as possible given the demands placed on information 
content. We make no claim that these data are com- 
plete in any volume limited sense. For our present 
purposes, this is less important than having data rep- 
resentative of all morphological types, luminosities, 
etc. This is because we arc more interested in the 
evolution of galaxies as physical systems than in the 
number density of any given type. Nonetheless, the 
potential effects of selection upon interpretation war- 
rant further discussion. 

The optical properties of the galaxies are plotted 
in Figs. 2 and 3. These illustrate some important 
basic facts about spiral galaxies. Disks exist over 
a large range in luminosity, surface brightness, and 
size. No particular value is characteristic for any of 
these; each has some distribution (McGaugh 1996a; 
de Jong 1996). Lower surface brightness disks tend to 
be bluer (Fig. 3), not redder as expected if they have 
fading evolutionary histories (McGaugh & Bothun et 
al. 1994). The lack of red low surface brightness galax- 
ies may itself be a selection effect (McGaugh et al. 
1996a), but one still wishes to understand why the 
ones we do know about are so very blue. 

Flux limited surveys are most efficient at identi- 
fying the brightest (L*) and highest contrast (/io*) 
galaxies that exist in the intrinsic distributions in a 
given passband. As a direct consequence, galaxy cat- 
alogs arc numerically dominated by such objects far 
out of proportion to their actual numbers. Samples 
selected for detailed study from such catalogs natu- 
rally tend to concentrate on these most prominent 
galaxies. We have avoided this by augmenting such 
data with our own data for low surface brightness 
galaxies. As a result, the combined data provide as 
comprehensive a view of spiral galaxies in the field as 
is currently possible, spanning the full range of disk 
morphologies from Sa to Sm/Im with /xq < 25 and 
h > 1 kpc. 

Though this procedure adds sorely needed weight 

to late types and low surface brightnesses, by no 
means does it guarantee a homogeneous sample. One 
indication of this is present in Fig. 2 as an apparent 
lack of low luminosity, high surface brightness galaxies 
with Mb > —18 and /zq < 23. This gap is not appar- 
ent when Ho is plotted against h (McGaugh 1996b), 
and does not necessarily indicate a real lack of such 
galaxies. In effect, the sample of de Jong contains 



only bright {Mb < ^18) galaxies while our own con- 
tains only ones of low surface brightness (/Uq > 23). 
The excluded region could well be populated, but no 
data which meet all our information requirements ex- 
ist in this region of parameter space. 

To test whether this can have a serious impact on 
our results, we perform the analysis described below 
for both the entire sample, and a subsample limited to 
Mb < — 18. Though the lever arm for fainter galax- 
ies is of course weaker, no serious differences occur 
among gas related parameters. The reason for this is 
simple: the galaxies are optically selected. Conclu- 
sions about the relation between optical parameters 
relevant to selection {Mb, Mo, and h) are difficult to 
draw in the absence of a rigorously complete sam- 
ple. However, the investigation into the evolution of 
galaxies is not much hindered because a wide range of 
these parameters are at least represented. Moreover, 
the initial selection of the galaxies is independent of 
Mhi- These data tell us quite a bit about the evo- 
lution of the galaxy types we do know about, and 
delimit what may still be missing. 

3. Gas Mass Fractions 

The baryonic gas mass fraction is 



where Mg is the total mass in the form of gas, and 
M* that in stars. To relate these to the observations, 
we need = and Mg = 77 Mhi- The parameter 
T* is the mass to light ratio of the stellar population, 
while T] plays an analogous role for the gas. It ac- 
counts both for helium and any mass in metals and 
for gas phases other than atomic. 

With these definitions, it is straightforward to ob- 
tain 

The flux ratio Mui/L is distance independent and has 

small errors in most cases. However, the conversion 
factors and rj must be estimated from other infor- 
mation. 

3.1. Stellar Population Mass to Light Ratios 

We wish to use the observable star light as an in- 
dicator of the stellar mass. Our stellar mass to light 
ratio T* includes all mass in stars living or dead since 
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we wish to distinguish between the gaseous and stel- 
lar mass reservoirs. The value of T* can be esti- 
mated from population synthesis models, but these 
are highly uncertain and depend on many factors like 
the IMF, star formation history, metallicity, gas recy- 
cling, the poorly known behavior of some late stages 
of stellar evolution, etc. Even worse for our purposes, 
it depends iteratively on the quantity of interest, fg. 
Though such models are still very useful, it is of obvi- 
ous interest to have independent means of constrain- 
ing T*. 

We use the dynamical constraints imposed by ob- 
servations of the vertical stellar velocity dispersion 
(Bottema 1993), as generalized by de Blok & Mc- 
Gaugh (1997). This gives 



T, = 1.936 X loO-'i(B-v) _ ;^_943 



(3) 



in the i?-band. This expression is derived from the 
disk Tully-Fisher relation with normalization imposed 
by the velocity dispersion observations. There is some 
uncertainty in the absolute normalization, but con- 
siderably less than in estimates lacking the boundary 
condition imposed by dynamical measurements. The 
color term has a residual dependence on population 
synthesis which goes in the expected sense that red- 
der galaxies have higher T*. This reproduces the ob- 
served trend well (Fig. 4) even though one does not 
expect B — y to be a perfect indicator of the variation 
of T* with the many underlying variables. 

The most important aspect of equation (3) is the 
normalization to dynamical constraints which give the 
most reliable estimate of the mass contained in stars. 
Without this, it is always possible to contrive models 
in which the stars contain all or none of the baryonic 
mass, and fg remains highly uncertain. With it, fg is 
reasonably well determined. 

The numbers returned by equation (3) are consis- 
tent both with stellar population models (e.g., van 
den Hoek et al. 1996) and rotation curve analyses 
(Sanders 1996; de Blok & McGaugh 1997) which are 
independent of the velocity dispersion measurements. 
As a further independent check of our results, we 
will also compute fg from a purely population syn- 
thesis estimate of in the /-band. The /-band 
is thought to be the most robust indicator of stellar 
mass (Worthey 1994), though one still does not ex- 
pect a universal value for all populations. The models 
of Bruzual & Chariot (1993) for a galaxy with con- 
stant star formation suggest a value of k, 1.2 over 
a wide range of ages, with fairly modest variation. 



Though both the absolute normalization and the uni- 
versal applicability of this value is uncertain, the use 
of this fixed number in the /-band gives indistinguish- 
able results from those of equation (3) in the B-band 
(Fig. 5). 

3.2. Gas Phases 

For the gas, 77 must be corrected for both the hy- 
drogen mass fraction X, and the phases of gas other 
than atomic. We assume a solar hydrogen mass frac- 
tion, giving T] = = 1.4. Though not universally 
correct, most of this is due to primordial helium. Vari- 
ations in helium and metal content result in devia- 
tions from this value of less than 10%. This is a very 
small effect compared to that of other gas phases. 
Ionized gas in H II regions and hotter plasma is of 
negligible mass in spiral galaxies, but not so molecu- 
lar gas. Hence, 



?7= 1.4 



1 + 



M(H2) 



M(HI) 



(4) 



The ratio of molecular to atomic gas varies systemat- 
ically with Hubble type (Young & Knezek 1989). We 
parameterize their mean relation by 



M(H2) 
M(HI) 



3.7 - 0.8T 0.043T2 



(5) 



(Fig. 6). For very late types, we have forced the fit 

to approach A/(H2)/M(HI) -> as T ^10. Equa- 
tion (5) gives a good fit to the trend of the data of 
Young & Knezek (1989) except perhaps for very early 
types, but there are few galaxies in the sample with 
T = 1. Of greater concern is the large amount of real 
scatter at each type — T is by no means a perfect 
indicator of M(H2)/M(HI), so equation (5) is only 
a statistical estimator of the relative molecular gas 
content. 

The normalization of equation (5) depends on the 
CO to H2 conversion factor, the precise value and 
universality of which is a matter of much debate. 
The conversion factor appears to be metallicity de- 
pendent (Wilson 1995), and later types are generally 
more metal poor than early types. If one were to cor- 
rect the noisy correlation of A/(H2)/M(HI) with T 
for the noisy correlation of Z with T, it would have 
only a modest effect on equation (5) since the inferred 
variation in the conversion factor with metallicity is 
fairly small compared to the range in Fig. 6. Little 
will change because early types are approximately so- 
lar metallicity and will receive no correction from the 
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canonical value of the conversion factor. The latest 
types are roughly 0.1 (McGaugh 1994a), corre- 
sponding to a factor of 5 increase in the conversion 
factor. However, M(H2)/M(HI) is inferred to be very 
small in late types, and a five fold increase in a small 
number still results in a very small correction to the 
final fg. Rather than convolve together all these un- 
certainties, we simply use equation (5), realizing that 
this is probably the greatest uncertainty in the pro- 
cess. 

The uncertainty in the correction for molecular gas 
will have little impact on our final results. The correc- 
tion is only important in early types, which are rela- 
tively gas poor fg ^ 0.3. It is the late types which are 
inferred to be gas rich, on the basis of the HI alone. If 
we are missing enormous amounts of molecular gas in 
late types it will only make yet more gas rich galax- 
ies which are already inferred to be very gas rich. 
Hence our qualitative results are very robust, even if 
the quantitative details are subject to improvement 
in the CO data. Moreover, it is only the ratio T*/r; 
that matters to our analysis (equation 2). This ra- 
tio does not vary much, and is not correlated with 
surface brightness. In order to affect our results, this 
ratio would have to vary systematically with surface 
brightness by an order of magnitude in just the right 
way to offset the strong observed trend in Mhi/L. 
Not only does this seem unlikely, it simply can not 
happen. The gas fractions of low surface brightness 
galaxies can not be lowered by large T, as these are 
not allowed by the rotation curve data (de Blok & 
McGaugh 1996), nor can ry be raised for high sur- 
face brightness galaxies where CO is readily detected. 
Hence the large observed variation in Mhi/L must 
correspond to one in fg. 

The role of molecular and atomic gas: It is 

generally thought that star formation occurs only in 
molecular gas. However, it appears that there is very 
little and perhaps no molecular gas in the lowest sur- 
face brightness galaxies (Schombert et al. 1990). We 
should therefore consider the possibility that stars 
form directly in from atomic gas. It is not obvious 
that this is necessary, though. The star formation 
rates are low, and a very modest amount of molecu- 
lar gas might exist as an intermediary step. Indeed, 
one could argue that the lack of molecular gas inhibits 
star formation and keeps the surface brightness low. 
While these are interesting possibilities, they do not 
matter to our present analysis, which is concerned 



only with the global supply of gas and not the details 
of its local consumption. In this context, it is per- 
haps best to think of the HI as a reservoir for future 
star formation, regardless of whether molecular gas is 
a mandatory intermediate step. Indeed, gradual ra- 
dial infall from the extended HI disks of spirals may 
be the source of the inferred accretion of unprocessed 
material since there is no evidence for substantial HI 
reservoirs in spherical halos (Bothun 1985). 

4. The Relations 

The gas mass fractions computed as described 
above are tabulated in Table 1 together with the in- 
put data. In general there is close correspondence 
between fg computed from S-band and /-band data. 
This is encouraging since independent methods are 
used to estimate which is the more important fac- 
tor in equation 2. Variation in the two different de- 
terminations of fg give some idea of the uncertainties 
therein. 

In Fig. 7 we plot the gas fractions against the B- 
band properties of the galaxies. A strong relation be- 
tween gas fraction and total luminosity, disk surface 
brightness, and color is apparent. The same relations 
with luminosity and surface brightness are also appar- 
ent in the /-band data (Fig. 8). There is no relation 
at all with the size of the optical disk. Fits to the 
i?-band data give 

fg = 0.12(Mb + 23) (6a) 
fg = 0.12(/io - 19) (6b) 
fg = -IA[{B -V)- 0.95] (6c) 

where the uncertainty in the slope is ±0.07 for Mb, 
±0.04 for /ig, and ±0.5 for color. These relations are 
of course subject to modest variation with improve- 
ments in the relations governing and r], but it is 
clear that dimmer galaxies are systematically more 
gas rich. The slopes of the fg-Ms and -/xq relations 
are indistinguishable. 

An important consequence of the systematic in- 
crease of gas fraction with surface brightness is the 
end of the gas consumption paradox. Since roughly 
equal numbers of galax;ies exist at each surface bright- 
ness, there are many galaxies which still retain sub- 
stantial gas reservoirs (see McGaugh 1996a for the 
distribution functions). Perhaps it is not surprising 
that the galaxies which are the most prominent op- 
tically are those which have converted most of their 
gas into stars. 
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4.1. Correlations 

Correlation coefficients from a principle component 
analysis of the S-band data are given in Table 2. 
When reading this table, be sure to remember the 
backwards convention in the definition of magnitudes. 
This affects the signs of the correlation coefficients: 
TZ = —0.59 between Mb and h means that brighter 
galaxies tend to be bigger, not smaller. 

First note that fg is well correlated with Mui/L, 
from which it is derived. This was not entirely a fore- 
gone conclusion, since the relations for T* and ry could 
conceivably have offset this. The use of these relations 
introduces noise which reduces the correlation coeffi- 
cient from unity. 

Amongst independent variables, the strongest cor- 
relation is between Mhi/L and /Uq. This relation is 
more significant than either of the well known rela- 
tions between Mhi/L and Mb (Rao & Briggs 1993; 
Salpeter & Hoffman 1996) or Mhi/L and T (Roberts 
1963), though these are also obvious in Table 2. This 
strong Mhi /L-^q relation translates into a strong one 
between fg and /xq. The central surface brightness of 
a galaxy is a good indicator of its cumulative gas con- 
sumption and hence its evolutionary state. 

The gas fraction is also correlated with Mb and 
B — V. Though both are used in the calculation of 
fg, this correlation is not an artifact. It stems from 
the basic fact that there is relatively more 21 cm lu- 
minosity in spirals of low optical luminosity, which 
also tend to be bluer. The relation with color goes in 
the sense that blue galaxies are more gas rich, as one 
would expect. Moreover, the same result is obtained 
when color is not used in the calculation of fg with 
the /-band data. 

The data appear to occupy a fairly narrow struc- 
ture in fg-HQ-color-MB space. In analogy to the 
Fundamental Plane for ellipticals, this may indicate 
the existence of an "Evolutionary Hypersurface" con- 
straining the evolution of spirals. Howevca\ this is 
violated by the extreme galaxy Malin 1 which is ex- 
ceedingly gas rich for its luminosity. A more funda- 
mental plane may be that in /g-/io-color space, but 
this too may be modified by further selection effects 
(see below). 

The gas fraction is not well correlated with size. 

The evolutionary rate of spiral galaxies is apparently 
scale free. Neither is fg correlated with Hubble type, 
in spite of the Mhi/L-T relation. This is because 
the dominant gas phase changes rapidly with Hubble 



type (Young & Knezek 1989), offsetting in the total 
gas fraction the trend apparent in the atomic phase 
alone. 

Among optical parameters, Hubble type depends 
most on Mb and /Uq. Though it is not surprising that 
these are relevant to the morphological appearance 
of a galaxy, T does not provide a clean measure of 
any physical parameter (Naim et al. 1995; McGaugh 
1995b; de Jong 1995). Type is not well correlated 
with size, though there is some tendency for earlier 
types to be redder. 

The linear size of a disk galaxy is not well corre- 
lated with anything except absolute magnitude. Ab- 
solute magnitude is appears to be reasonably well cor- 
related with surface brightness in the obvious sense 
that higher surface brightness galaxies are brighter. 
However, this may be illusory, being the relation most 
subject to selection effects. In particular, essentially 
all the low luminosity galaxies are low surface bright- 
ness (recall Fig. 2). This may or may not be represen- 
tative of reality. If we restrict the sample to objects 
brighter than —18 where the Mb-Ho plane is more 
completely covered, this correlation coefficient drops 
to TZ = 0.34. That between /xq and h increases to 
TZ = 0.44 — bigger galaxies tend to be lower surface 
brightness. These are the only significant changes in 
the correlation matrix when the sample is restricted 
to bright galaxies. 

Optical galaxy selection depends jointly on Mb, 
Ho, and h (Disney & Phillipps 1983; McGaugh et 
al. 1995a), so correlations between these variables are 
very sensitive to sample selection and must be treated 
with caution. Correlations with other variables are 
not as directly subject to selection effects, and can 
sensibly be analyzed as long as representative data is 
available over a reasonable dynamic range. In partic- 
ular, correlations involving fg appear to be robust. 

4.2. Further Selection Effects 

We should nevertheless consider whether the rela- 
tions derived above are affected by selection effects 
which in the past helped conceal them. We have ad- 
dressed the problem of isophotal selection, to the ex- 
tent currently possible, by specifically targeting low 
surface brightness galaxies. This increases the dy- 
namic range available in /kq, but it is obvious that 
such effects remain from the truncation of the data 
as /io 25. The galaxies studied here are selected 
optically, but of course a galaxy must have detected 
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HI to perform this analywis. The lack of low gas frac- 
tion, low surface brightness galaxies could therefore 
be a further selection effect. 

The relations appear continuous across samples, 
which would seem to argue against strong selection 
effects in this sense. The 21 cm detection rate is very 
high (80 %) for low surface brightness galaxies in the 
Schombert et al. (1992) catalog, so we do not seem 
to be missing a significant number of these predom- 
inantly spiral galaxies. However, the fact that these 
galaxies are very blue (Fig. 3) may itself be a selec- 
tion effect, since they were selected on blue sensitive 
plates (see McGaugh 1996a). For this to be impor- 
tant to Fig. 7, a large additional population of gas 
poor red low surface brightness galaxies would also 
have to exist. 

Other types of galaxies (e.g., dE and dSph) are cer- 
tainly known to exist which have low surface bright- 
nesses and low gas fractions. These particular types 
bear no apparent relation to the larger low sur- 
face brightness spirals cataloged by Schombert et al. 
(1992), or to higher surface brightness spirals. Never- 
theless, it may well be that the observed relations rep- 
resent only upper envelopes, and a substantial popu- 
lation of gas poor low surface brightness galaxies re- 
mains to be discovered. Indeed, such a popiilation 
may be demanded by the rapid evolution of the lu- 
minosity density (Dalcanton 1993; Lilly et al. 1995 — 
see §6). The lack of gas rich, high surface brightness 
galaxies is not a selection effect: such objects would 
be prominent both optically and in HI. 

5. Disk Evolution 

Irrespective of these caveats, clear evolutionary 
trends exist. The lower surface brightness and lower 
luminosity spirals that wc know about have consumed 
less of their gas and are less evolved than brighter sys- 
tems. This result is consistent with other lines of evi- 
dence from star formation thresholds (van dcr Hulst et 
al. 1993), abundances (McGaugh 1994a; Ronnback 
& Bergvall 1995), colors (McGaugh & Bothun 1994; 
Ronnback & Bergvall 1994; Bergvall & Ronnback 
1995; de Blok et al. 1995) and environments (Bothun 
et al. 1993; Mo et al. 1994), all of which indicate that 
low surface brightness galaxies evolve slowly, and per- 
haps form late. 

On their own. Figs. 7 and 8 suggest a straight- 
forward picture: disk galaxies form and then evolve 
at constant size, converting gas into stars and hence 



changing their luminosity and surface brightness at 
the same rate. This assures strong fg-Ms and fg-j-io 
relations with the same slope. There need be no cor- 
relation between Mb and /io since galaxies exist over 
a range of sizes. Disk galaxy evolution thus seems 
to proceed in an orderly fashion, without excessive 
amounts of merging, size evolution, strong stochastic 
star bursts, or other exotic phenomena. 

The evolutionary trajectory of a galaxy in Fig. 7 
or 8 depends on the evolution of stars it contains, and 
the rate at which those stars are made. The star for- 
mation rate (SFR) is just the negative of the rate of 
gas consumption: = —fg. at least mitil gas recy- 
cling becomes important at low fg. The photometric 
evolution of stars depends on many factors (e.g., their 
metallicity) which lead to a variety of complications, 
but the basic physics is well understood and allows 
one to estimate the luminosity that will be produced 
by a given mass of stars. This gives the abscissa. 
The ordinate follows from the SFR which is tradi- 
tionally parameterized as some convenient function 
of time (e.g., Guideroni & Rocca-Volmerange 1987). 
For declining SFRs in which most of the star forma- 
tion occurred well in the past, a commonly assumed 
form is 

f^ = -fg = r-'e-'^\ (7) 

The SFR is normalized to the total baryonic mass of 
a given galaxy so that < /, < 1. The e-folding 
time scale r is typically several Gyr as thought ap- 
propriate for early type galaxies. For SFRs which are 
constant (appropriate for late type spirals) or rising 
(as is perhaps the case for irregular galaxies), a useful 
form is 

In this case, Tg is the time to total gas consumption. 
This is not quite the same as an e-folding time scale, 
but in both cases sets the evolutionary rate. 

We do not delude ourselves that either of these sim- 
ple SFRs accurately represent what really occurs in 
any given galaxy. They do nevertheless provide a very 
useful way of parameterizing the range; of possible av- 
erage star formation histories. This simple analytic 
approach does enable one to infer several interesting 
constraints on the more complicated reality when a 
large dynamic range in fg is available. 

Integration of these SFRs allows us to define sev- 
eral quantities of interest (and also illustrates why 
these forms are simple and popular). The age of a 
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galaxy. Tg- wc; will take to mean the time; since star 
formation commenced. This is the obvious working 
definition of galaxy formation in this context. It does 
not. of course, specify precisely how a galaxy formed. 
Star formation might well begin while the collapse of 
gas is ongoing, or might have begun in separate sub 
units which subsequently merge into the single entity 
observed today. It is convenient to imagine the case 
in which the time scale for such formation events is 
short compared to the evolutionary time scale r so 
that formation is contemporaneous with the onset of 
star formation and Tq is the true age of the galaxy 
for all practical purposes. This need not be the case, 
so one should keep in mind that Tq really measures 
the time since the onset of star formation. 

Another interesting quantity is the nic;an age of the 
stars, (T*). It is often possible to vary the shape of the 
SFR so as to arrive at any desired fg with a degener- 
ate set of ages Tq and evolutionary rates r. However, 
otherwise plausible combinations of these two some- 
times have implausible consequences for (T*). These 
and other quantities and limits of interest are given 
in Table (3). 

Equipped with this formalism, we are nearly ready 
to interpret the positions of galaxies in Fig. 8. The 

value of fg comes from stipulating an SFR, but the 
precise trajectory a galaxy follows in luminosity de- 
pends on the details of the stars for which a popula- 
tion model is required. This gives the variation of T* 
with time. The amount of light per unit mass pro- 
duced by a single generation of stars declines rapidly 
as luminous, massive stars evolve and die. However, 
if star formation is ongoing for any period of time, 
this effect is offset by the birth of new stars. To esti- 
mate the interpay of these effects, we need to choose 
a specific population model. For now, we employ 
the models of Guidcroni & Rocca-Volmerange (1987), 
briefly described in Table (4). This particular choice 
is made for the simple reason that they published the 
full range of fg for an interesting variety of SFRs. Var- 
ious improvements in modeling have since been made, 
but no uniform set of models is readily available over 
the entire dynamic range we need. Moreover, we are 
trying here to avoid as much as possible the serious 
astronomical problem of model dependency; specific 
models are adopted only to illustrate the various is- 
sues. Compared to the dynamic range in fg, stellar 
populations experience fairly mild evolution of in 
all plausible models. So while the detailed shape of 
evolutionary tracks is highly model dependent, the 



basic gist is not. 

For example. Fig. 8 by itself does not represent a 
direct evolutionary sequence. Galaxies do not start 
at faint M and /xq and simply brighten along the 
observed sequence. The slope of the observed rela- 
tions is an order of magnitude different from what 
is expected from evolutionary models for the rate of 
change of absolute magnitude (Am) with gas fraction 
(A/g). Instead of the slope observed in equation (6), 
= 8, a very broad range of models evolve slowly 
in luminosity compared to the rate of gas consump- 
tion with -2 < -f^ < 2 (Fig. 9). This difference 
between observed and model slopes requires a sys- 
tematic variation in the evolutionary rate with sur- 
face brightness. Progressively dimmer galaxies either 
evolve more slowly or are younger, or some combina- 
tion thereof. 

The isochrones of the models of Guideroni & Rocca- 
Volmerange (1987) are roughly parallel to the ob- 
served slope for ages Tg > 10 Gyr. This is 
achieved by systematically changing the gas consump- 
tion time scale with Hubble type (Table 4). Since 
very little of these data were available when the mod- 
els were constructed, this is a success of the standard 
picture of galaxy evolution the models represent. On 
the other hand, the correlation between M-m/L and 
type was known, and the models were constructed to 
match that data. 

From Fig. 9 one can read off the combined age- 
SFR of any given galaxy. If one can be constrained by 
other information like colors, the other follows. Type 
is correlated with surface brightness (Table 2; de Blok 
et al. 1995; de Jong 1996), so it is possible to make 
all spirals the same old (~ 15 Gyr) age and satisfy 
the observed gas fractions simply by varying the time 
scale T. The variation must be quite systematic with 
surface brightness, and so large that one has to change 
the functional form of the SFR to accommodate it. 
(In principle, one could retain an exponential SFR by 
allowing t ^ oo and then becoming negative.) 

The spread of the data is large compared to the 

separation of the isochrones. This may indicate a 
large dispersion in age. More likely, it simply results 
from an overlap of the range of baryonic masses of 
any given SFR type. That is, a particular SFR type 
like Sd may be experienced by galaxies of a range of 
masses, not just one particular one as drawn in Fig. 9. 
Because there is inevitably overlap of this sort we do 
not attribute any special significance to the kink in 
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the model grid for Sa and Sb SFRs at low fg — the 
overlap would wash it out, even if the models were 
credible at this level of detail. The lack of any cor- 
responding feature in the data might indicate that t 
does not vary as widely as assumed in the models, so 
that for spiral disks deviations from a constant SFR 
are not enormous (r ^ 8 Gyr). 

The standard picture of galaxy evolution repre- 
sented by the models of Guideroni & Rocca-Volmerange 
(1987) must attribute the observed change of fg with 
surface brightness to variations in t because Tq is 
taken to be the same for all galaxies. Once this ar- 
tificial boundary condition is imposed, there is no 
choice. It is not obvious that this is really a sen- 
sible assumption. If galaxy formation occurs from 
the gravitational collapse of primordial fluctuations 
in the density field, then the epoch of collapse of any 
given object depends on the density contrast of the 
fluctuation from which it arises. Most of the observa- 
tional evidence (McGaugh 1992; 1996b) suggests that 
this occurs in the obvious way: large density con- 
trasts result in high surface brightness galaxies, and 
small density contrasts result in low surface bright- 
ness galaxies. If this is indeed the case, one expects 
lower surface brightness galaxies to be younger as in- 
ferred from their colors (McGaugh & Bothun 1994; de 
Blok et al. 1995). The absolute magnitude of the age 
difference may or may not be large; it is very difficult 
to put an absolute number on. From the colors, the 
relative shift typically appears to be on the order of 
a few (3 - 5) Gyr. If high surface brightness spirals 
are 13 or 14 Gyr old, low surface brightness disks are 
roughly 8 or 10 Gyr old. 

To examine this possibility in the context of the 
gas fraction data, consider as a trivial limiting case a 
scenario in which all spiral disks have the same shape 
star formation history. The evolutionary trajectories 
would all look the same in Fig. 9, and (by choice) 
isochrones would be horizontal lines. Instead of a sys- 
tematic variation in the evolutionary time scale, real 
age differences woiild be required to explain the trcmd 
of the data. For the ease of constant star formation at 
the same relative rate (the Sd track), bright galaxies 
with fg Ri 0.3 would have ages Tq w 17 Gyr, while 
the dim ones with fg « 0.6 would be Tq « 10 Gyr 
old. 

Either of these scenarios seem plausible, and real- 
ity is presumably more complicated — both age and 
evolutionary rate may play a role. Gas fractions alone 
do not allow us to distinguish which factor dominates. 



They do show that clear differences exist. This may 
help lift the degeneracy that has plagued analyses of 
colors alone. If so, it may become possible to age-date 
spiral galaxies. 

The distinction between a galaxy which is gas rich 
today because it formed early and evolved slowly or 
because it formed late is not clear. Some limit can 
be placed on these extreme cases by the plausibil- 
ity of the required SFR and mean stellar age. If we 
insist that all galaxies have the same age, the only 
way to have very gas rich blue galaxies {fg > 0.4) at 
= is to have star formation rates which increase 
strongly with time. These retain much gas today sim- 
ply because very little happened for the first half of 
a Hubble time. For the Im track in Fig. 9, only 10% 
of the gas is consumed in the first 9 Gyr. The dis- 
tinction between this and a galaxy which forms late 
is inobvious, and returns us to the thorny issue of the 
definition of galaxy formation. Perhaps a lengthy col- 
lapse by gradual dissipation and a slow evolution are 
intimately related. 

For a more quantitative illustration, consider the 
conservative case of a declining star formation rate. 
Retaining a gas fraction greater than what is con- 
sumed in one e-folding {fg > = 0.37, correspond- 
ing to fiQ > 22) requires an evolutionary time scale 
T > Tu greater than the age of the universe, with a 
mean stellar age (T,) < O.GTu. This follows simply 
from requiring that the galaxy be younger than the 
universe: To < Tu- If we adopt a rising star forma- 
tion history as suggested by the Im model in Fig. 9, 
the evolutionary time scale can be very (but not ar- 
bitrarily) long. This is only achieved at the expense 
of making the stars very young: (T*) < ^/iTu, i.e., 
(T*) = 3 Gyr for Tq = 12 Gyr. Either way, most of 
the star formation in very gas rich galaxies must be 
weighted towards more recent epochs (cf. McGaugh 
& Bothun 1994); irrespective of whether the galaxy 
formed early and did nothing, or formed late. Indeed, 
galaxies with fg > 0.5 have yet to realize most of their 
potential for star formation, which must be weighted 
towards future epochs. 

6. Conundrum 

The picture which emerges thus far seems fairly 
sensible, confirming and extending established wis- 
dom about the evolution of local galaxies. Related in- 
formation can be inferred from deep redshift surveys, 
and considerable progress has been made in measur- 
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ing the properties of galaxies over an interesting range 
of redshift. These studies indicate a large amount of 
evolution in both field (Lilly et al. 1995; Ellis et al. 
1996) and cluster (Schade et al. 1996) galaxies. High 
resolution imaging with HST shows that most of this 
evolution is due to morphologically late type galaxies 
(Glazebrook et al. 1995; Driver et al. 1995, Abraham 
et al. 1996). This is consistent with the evolution of 
the luminosity function constructed from deep sur- 
veys: red galaxies appear to evolve little, while blue 
ones have been fading very rapidly (Lilly et al. 1995; 
Ellis et al. 1996). 

There is good evidence for direct evolution of the 
surface brightnesses of bright disks (Schade et al. 
1995, 1996). These appear to have surface bright- 
nesses which are ^ 1 mag arcsec"^ brighter at z « 0.5 
than at z = 0. This is perfectly consistent with 
Fig. 9: these galaxies were presumably more gas rich 
at z = 0.5 and have since evolved along a track like 
that plotted for early-type (Sa-Sc) spirals. 

The local and high redshift approaches to the evo- 
lution of early type spirals therefore give nicely con- 
sistent results. The same can not be said for late 
types. Many studies now indicate that the excess faint 
blue galaxies are predominantly late types which are 
evolving rapidly. If one looks locally for faint, blue, 
late type galaxies, one finds low surface brightness 
galaxies. Indeed, the similarity between the physical 
properties of these two populations led to the sugges- 
tion that one might account for the other (McGaugh 
1994b). This similarity now extends to morphology: 
these galaxies literally look alike (cf. Dalcanton & 
Shectman 1996). However, a simple one to one con- 
nection is not sufficient: although deep surveys are 
more sensitive to low surface brightness galaxies than 
local ones, they are still selected by flux. Magnitude 
limited surveys are always numerically dominated by 
the highest surface brightness galaxies regardless of 
how common low surface brightness galaxies actually 
are (McGaugh et al. 1995a), so for most plausible sce- 
narios low surface brightness galaxies do not make a 
significant contribution to the faint counts (Ferguson 
& McGaugh 1995; McLeod & Rieke 1995) unless they 
evolve siibstantially. 

It is tempting to conclude that this is exactly what 
happens: late type blue galaxies evolve rapidly, fad- 
ing to become low surface brightness by the present 
epoch (Phillipps & Driver 1995). While sensible in 
the context of the high redshift data, it makes no 
sense for the known local galaxies. Fading scenarios 



for the non-dwarf low surface brightness galaxies in 
the present sample are strongly excluded by both the 
colors (McGaugh & Bothun 1994; Fig. 3) and by the 
gas fractions. In order to have a strongly fading tra- 
jectory (Sa or even more extreme) in Fig. 9 and retain 
the high observed gas fractions, such galaxies would 
have to be extremely young: Tq < 2 Gyr for r = 3 
Gyr and fg > 0.5. They would not have even formed 
yet at the redshifts where the inferred evolution oc- 
curs. 

Indeed, the evolutionary histories discussed in §5 
are precisely the opposite: late type, low surface 
brightness galaxies evolve slowly, not rapidly. Photo- 
metrically, they evolve at roughly constant luminosity 
or even brighten with time. We are left with a puz- 
zling c;ornindrum: there exist two populations, the 
faint blue galaxies at high redshift and low surface 
brightness galaxies locally, which have very similar 
colors, clustering properties, and morphologies: they 
literally look alike, yet they can not be the same. 

Perhaps the picture which emerges from the deep 
survey work is deceptively simple. Dividing the lumi- 
nosity fimction into bins by color (Lilly et al. 1995) or 
spectral type (Heyl et al. 1996) gives the appearance 
of a slowly or even non-evolving red population and 
a rapidly evolving blue one. Yet one expects rather 
more evolutionary fading for red galaxies, consistent 
with the amount of surface brightness evolution ob- 
served by Schade et al. (1996). Perhaps there is some 
cross-over between color bins. This is expected at 
some level, as even early type galaxies should have 
been blue at some point early in their evolution. 

It is not obvious that cross-over is adequate to ex- 
plain the evolution of early types, but it may well 
be (Pozzctti et al. 1996). A more serious problem 
is understanding the rapid evolution of the blue lu- 
minosity density attributed to late types. Since the 
faint blue galaxies we do know about locally seem un- 
able to explain this, perhaps the inference of a rapidly 
evolving dwarf population is correct and there exists 
locally an additional remnant population (Phillipps & 
Driver 1995; Babul & Ferguson 1996). Such a popu- 
lation would presumably now have faded and become 
red, low surface brightness, and gas poor. This is pre- 
cisely the part of the fg-fio plane not constrained by 
observation. A very large population of such objects 
woiild be required, but it is quite possible that such 
a population would have remained undetected so far 
(McGaugh 1996a), especially if the typical object is 
intrinsically small as indicated by deep HST imaging 
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(Im (it al. 1995). This picture might bo considered at- 
tractive for avoiding the reverse of the gas consump- 
tion paradox: if gas rich galaxies exist as well as those 
near to depleting their gas supply, why arc there no 
galaxies well past this point? The unavoidable conse- 
quence of such a scenario is that our knowledge of the 
local galaxy population remains glaringly incomplete. 

7. Conclusions 

Examination of the gas mass fractions of spiral 
galaxies over a large range in luminosity and surface 
brightness provides a number of insights into galaxy 
evolution. 

1 . The fraction of baryonic mass in spiral galaxies 
which has been converted from gas into stars is 
correlated with luminosity, surface brightness, 
and color. It is not correlated with the scale 
length of the disk. 

2. The strongest correlation is between gas con- 
tent and the central surface brightness of the 
disk. The surface brightness of a disk is a fun- 
damental parameter and a good indicator of the 
evolutionary status of a galaxy. Dim galaxies 
evolve slowly, have relatively young stellar pop- 
ulations, and may have formed late. 

3. Since luminosity traces baryonic mass and sur- 
face brightness traces sm-face mass density, the 
strength of the /g-/io relation STiggests that s\ir- 
face density is at least as important as total 
mass in determining the evolution of a galaxy. 

4. The fg-Ms and fg-Ho relations have the same 
slope. Any size evolution would cause surface 
brightness evolution over and above that due to 
luminosity evolution, leading to a difference in 

these slopes. The lack of such a difference indi- 
cates that spirals evolve at roughly fixed size. 

5. The large number of gas rich low surface bright- 
ness galaxies resolves the gas consumption para- 
dox. 

6. The populations of local low surface brightness 
galaxies and high redshift faint blue galaxies 
are physically and morphologically similar yet 
grossly different in the way they evolve. This 
poses a difficult conundrum. One possible so- 
lution is an additional population of local, red, 
gas poor low surface brightness galaxies. 



There is no reason to believe that we have reached 
the faintest limits of surface brightness, or that yet 
more gas rich galaxies do not exist. Given the bar- 
rier imposed by isophotal selection limits, there could 
certainly be more very low surface brightness galaxies 
that have yet to be discovered. Some of these could be 
even more gas rich (by extrapolation of the observed 
trends) and others might be gas poor (§6). 

Several further lines of work are suggested. 

1 . Further surveys for low surface brightness galax- 
ies, both to deeper isophotal limits and in red- 
der passbands. Blind 21 cm surveys with optical 
follow-up to measure surface brightnesses would 
also be interesting. 

2. Detailed modeling of individual galaxies to match 
colors, metallicities, star formation rates and 
gas mass fractions. 

3. Examining these properties on a local as well as 
global basis. 

These projects would address some of the puzzles 
posed, and further advance our understanding of 
galaxy evolution. 

In essence, we must understand how low surface 
brightness galaxies can remain gas rich until the 
present epoch. Either quiescence or youth, or some 
combination thereof, can be invoked. If quiescence is 
primarily responsible, it must be possible to suppress 
large scale star formation in gas disks for many Gyr, 
perhaps through critical density dependent phenom- 
ena (e.g., Kennicutt 1989). If youth is a factor, the 
existence of such unevolved galaxies suggests that the 
process of galax;y formation is ongoing, or has ended 
only recently. 

We are grateful to J. M. van der Hulst and L.B. 
van den Hoek for helpful discussions, and thank the 
referee for suggesting a more expanded discussion. 
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Table 1 — Continued 



Galaxy 


Mr 




h 




B - 


- V 


r 


Mhi 
Lb 


Jg 


Mhi 
Li 


Jg 


B 


- 1 


UGC 10584 


-21.20 


21.77 


8, 


.3 


0, 


.61 


5 


0.55 


0.47 


0.38 


0.44 


1, 


.76 


UGC 11628 


-21.36 


22.27 


11 


.0 


0, 


.96 


2 


0.26 


0.29 


0.18 


0.36 


1, 


.76 


UGC 11708 


-20.65 


21.51 


5, 


.0 


0, 


.66 


5 


0.27 


0.25 


0.16 


0.21 


1, 


.87 


UGC 11868 


-18.56 


22.65 


1, 


.7 


0, 


.74 


10 


0.42 


0.22 










UGC 11872 


-19.46 


20.46 


1, 


.5 


0, 


.84 


3 


0.21 


0.24 


0.11 


0.23 


2, 


.03 


UGC 12151 


-17.76 


23.27 


3 


.1 


0. 


.55 


10 


1.45 


0.59 


1.47 


0.61 


1. 


.33 


UGC 12343 


-20.93 


21.95 


8, 


.3 


0, 


.98 


5 


0.26 


0.18 


0.14 


0.22 


1, 


.98 


UGC 12379 


-21.37 


21.99 


8, 


.0 


0, 


.80 


4 


0.42 


0.36 


0.20 


0.31 


2, 


.14 
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-20.47 


21.50 


5. 


.1 


0. 


.64 


5 


0.73 


0.52 


0.48 


0.48 


1. 


.79 


UGC 12511 


-19.76 


22.46 


5, 


.8 


0, 


.72 


6 


1.52 


0.61 


0.76 


0.55 


2, 


.09 


UGC 12614 


-20.91 


20.94 


5, 


,1 


0, 


.73 


5 


0.20 


0.20 


0.14 


0.22 


1, 


.71 


UGC 12638 


-20.72 


22.17 


9, 


.1 


0, 


.75 


5 


0.30 


0.25 


0.21 


0.26 


1, 


.76 
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-20.17 


21.76 


5, 


.3 





.78 


4 


0.44 


0.38 


0.31 


0.43 


1, 


.72 
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-18.92 


24.59 


8, 


,4 


0, 


.37 


9 


1.28 


0.67 


1.31 


0.60 


1, 


.32 
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-18.74 


21.81 


3, 


,1 


0. 


.54 


6 


0.37 


0.35 


0.36 


0.36 


1, 


.37 


UGC 12776 


-21.25 


23.39 


20 


.3 






3 


0.38 




0.29 


0.45 


1, 


.62 


UGC 12808 


-21.41 


20.31 


3 


.9 


0, 


.70 


3 


0.09 


0.16 


0.05 


0.15 


1, 


.85 


UGC 12845 


-20.19 


22.69 


8. 


.0 


0. 


.89 


7 


0.79 


0.33 


0.56 


0.42 


1. 


.72 



Table 2 

Correlation Matrix 





/. 


Mb 


Mo 


log(/i) 


B-V 


T 




0.86 


0.58 


0.73 


0.00 


-0.37 


0.64 


fg 
Mb 

Mo 

log(/i) 




0.51 


0.63 
0.58 


0.03 
-0.59 
0.18 


-0.55 
-0.49 
-0.35 
0.33 


0.36 
0.68 
0.59 
-0.37 



B-V -0.45 
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Table 3 
Star Formation Histories 



Exponential^ Power Law'' 



Star Formation Rate'^''^ /„ = -r-^e"*/^ -2=±l(l.Y 

Gas Mass Fraction fg= e"*/"^ l-^^:^)"" 

Galaxy Age^ Tg = -rlnfg Tg(l - /g)V(^+i) 

Mean Stellar Age (T*) = -r (l + ^) Tg/{x + 2) 

Time Scale Limit<='f T>-Tu/l-afg < T,7(l - /g)-i/(^+i) 

Mean Stellar Age Limit^ (T*) < Tu (jj^ + Tu/{x + 2) 



''Declining SFR with r > 0. 

'Tlat or rising SFR with Tg and x > 0. 

'^T = e-folding time scale; Tg = gas consumption time scale. 

''Global SFR normalized by total baryonic mass. 

''Tg defined by the onset of star formation = age of oldest stars. 

^Limits obtained by requiring Tg <Tu, the age of the universe. 



Table 4 
Population Models 



T 


r (Gyr) 


X 


Form of SFR 


Sa 


3.3 




Exponential decay 


Sb 


5 




Exponetial decay 


Sc 


10 




Exponetial decay 


Sd 


20 





Constant 


Im 


20 


2 


Rising power law 



Note. — Hubble type is used here to de- 
note the models of Guideroni & Rocca- 
Volmerange (1987). Strictly speaking, these 
only specify the shape of the model star for- 
mation history. 
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Fig. 1. — The HI mass to light ratio as a function 
of disk central surface brightness, a) i?-band b) /- 
band. There is a strong correlation between the gas 
content of spiral galaxies and their surface brightness. 
This is apparent in both passbands even though there 
are fewer data in /. What scatter there is reflects 
real differences between galaxies as the uncertainties 
are small compared to the range plotted: typically 
0.2 - 0.3 mag. in /zo, and 0.2 dex in the distance 
independent flux ratio My{i/Lb- Note the sharp end 
to the data as iiq ^ 25 B mag arcsec"^ imposed by 
isophotal selection limits. Without our data for low 
surface brightness galaxies, this occurs brighter than 
1^0 ~ 23 in (a). This is also apparent in the gap in 
(b) at /Lto^ ~ 21. Examining only galaxies brighter 
than this conceals the correlation for lack of dynamic 
range. The very low surface brightness giant Malin 
1 [the triangle in (a)] breaks this isophotal barrier 
because it was scrcndipitously discovered, in part due 
to its high gas content. 

Fig. 2. — The optical (S-band) properties of spiral 
galaxies: a) central surface brightness vs. absolute 
magnitude and b) the scale length length vs. abso- 
lute magnitude. Different symbols represent different 
Hubble types, progressing from early to late types in 
order of increasing number of points on the symbol 
[key inset in (a)]. 

Fig. 3. — /-band central surface brightness /io^ = 
jiQ — {B — I) vs. V — I color. Symbols correspond to 
Hubble type as per Fig. 2. Note that lower surface 
brightness disks are systematically bluer than those 
of higher surface brightness — they are not the result 
of fading after cessation of star formation. 

Fig. 4. — The vertical velocity dispersion of spiral 
disks function of color. Data are from Bottema 
(1993; 1995). The line is our estimator for T* mapped 
onto the velocity dispersion by equation (8) of Bot- 
tema (1993). The normalization is chosen to match 
the data, but the shape of the line is not a fit thereto 
(see de Blok & McGaugh 1997). Nevertheless, it does 
match the observed trend, which goes in the expected 
sense that redder galaxies have higher T*. 

Fig. 5. The stellar mass of sample galaxies (in so- 
lar masses) computed from = T*L independently 
in the B (abscissa) and / (ordinate) bands. Sepa- 
rate measures of the luminosity {Lb and Lj) are used 
together with the estimators for discussed in the 
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text. The data closely follow the line of unity, indicat- 
ing reasonable agreement between the two methods. 

Fig. 6. — The ratio of molecular to atomic gas mass 
as a function of Hubble type. Data (mean and error 
in the mean) are from Young & Knezek (1989). As 
usual, the error in the mean understates the amount 
of real scatter at a given type. The line is our fit to 
the data (see text). 

Fig. 7. The gas fraction of spiral galaxies plotted 
against a) absolute B magnitude, b) disk central sur- 
face brightness, c) disk scale length and d) color. The 
full possible range of < /g < 1 is plotted and for 
the first time is needed — the dimmest galaxies are 
mostly gaseous with /g ^ 0.5. A strong relation exists 
between fg and Mb, j-J-a, and color. There is no rela- 
tion at all between fg and the size of the disk. Note 
that the fg-Ms and -/xo relations have the same slope. 
Malin 1 (triangle) is an extreme case with a gas frac- 
tion which is abnormally high for its luminosity but 
normal for its surface brightness. 

Fig. 8. — Like Fig. 7 except for the 7-band data. The 
figure shows fg against a) absolute / magnitude, b) 
/-band disk central surface brightness, c) scale length, 
and d)V — I color. The gas fraction is determined in- 
dependently from that in Fig. 7 by assuming that T{ 
is constant for all spirals. This gives indistinguishable 
results from the B-band even though a completely 
different estimator for the stellar mass has been em- 
ployed. 

Fig. 9. — The data from Fig. 7 (a) and (b) ovcrplotted 
after shifting by the intercepts of equation (3). Cir- 
cles: Mb- Squares: /kq. Lines are model evolutionary 
tracks from Guideroni & Rocca-Volmerange (1987). 
Horizontal lines are isochrones every 2 Gyr from 1 to 
15 Gyr. Vertical lines are evolutionary trajectories for 
specific star formation histories, labeled by the Hub- 
ble type the model attempts to represent. Early types 
(Sa to Sc) have exponentially declining star formation 
rates, while Sd is constant and Im rising with time 
(Table 4). The position of these lines depends on the 
total baryonic mass of a galaxy. They are drawn with 
separations corresponding to the mean of each type as 
stipulated by Guideroni & Rocca-Volmcrangc (1987); 
the absolute position of the resulting grid has been 
adjusted to coincide with the data. Though these (or 
any) population synthesis models are by no means 



perfect, tlicjy do illustrate how fg constrains the com- 
bination of evolutionary rate and age. With a large 
dynamic range in fg and the traditional constraints of 
color and metallicity, it may be possible to construct 
a multidimensional extragalactic analog to the HR 
diagram which might enable the age-dating of spiral 
galaxies. 
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